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ABSTRACT 


The  experimental  set-up  for  laboratory  study  of  spacecraft  control/structural 
interaction  has  been  designed.  Design  specifications  have  been  derived,  and  all 
the  actuators  and  sensors  have  been  selected  except  the  end-point  displacement 
sensing  of  the  arm.  The  mainbody  and  the  flexible  arm  have  been  fabricated  to 
meet  design  criteria. 

The  equations  of  motion  for  the  experimental  model  have  been  derived 
and  natural  frequencies  determined.  The  natural  frequencies  of  the  fle.xible  arm 
have  been  determined  experimentally  and  compared  with  analytical  predictions 
obtained  by  using  the  GIFTS  finite  element  analysis  program.  The  experimental 
and  analytical  results  are  in  good  agreement  except  the  first  mode. 
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I.  INTRODUCTION 


A.  THE  CONTROL  OF  FLEXIBLE  STRUCTURES 

Flexible  structures  have  become  an  integral  part  of 
modern  spacecraft  design  for  a  variety  of  reasons. 

Deployable  structures  extending  several  meters  into  space 
used  for  solar  arrays  are  required  for  todays  high  powered 
satellites.  Radar  reflectors  of  large  diameter  are  used  for 
communication  satellites.  Meanwhile,  as  the  systems  expand 
in  size,  weight  constraints  continue  to  plague  the  satellite 
designer.  Hence,  lighter  and  flexible  structures  are  often 
used.  This  results  in  lower  structural  natural  frequencies. 

The  large  antenna  reflectors  require  higher  attitude 
pointing  accuracy.  In  order  to  provide  higher  attitude 
pointing  accuracy,  the  control  bandwidth  has  to  be 
increased.  Because  of  the  decrease  in  the  structural 
natural  frequencies  and  increase  in  control  bandwidth,  there 
is  a  high  possibility  of  the  structural  frequency  falling 
into  the  control  bandwidth,  resulting  in  control/structure 
interaction. 

For  some  space  applications,  the  flexible  structure  will 
require  active  control.  As  a  result,  current  and  future 
spacecraft  study  of  control/structure  interaction  has  become 
a  challenging  problem. 
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B.  BACKGROUND. 


The  control  of  rigid  body  dynamic  systems  has  been  a 
subject  of  study  for  many  years.  Controlling  three-axis 
stabilized  satellites  has  been  successfully  achieved  for 
several  years.  The  1990 's,  however,  bring  a  new  challenge 
to  the  control  problem.  Future  NASA,  military,  and 
commercial  space  missions  will  involve  advanced  space 
systems  which  have  higher  power  reguirements ,  greater 
reguired  pointing  accuracies,  and  faster  slewing  through 
larger  angular  motions.  Additionally,  permanent  space 
stations,  like  Freedom,  and  complex  space-platforms,  like 
the  Hubble  Telescope,  will  reguire  the  use  of  large, 
flexible  manipulators  for  maintenance,  on-orbit 
construction,  etc. 

Research  into  the  demanding  problems  of  the  control  of 
flexible  structures  has  been  undertaken  for  the  last  two 
decades  by  major  space-system  engineering  firms,  the 
government  laboratories,  and  universities.  The  effort  has 
been,  however,  limited  to  analytical  studies.  Recently, 
experimental  work  has  started  at  some  universities. 

Texas  A&M  University  has  constructed  a  hub-appendage 
configuration  to  perform  large  angle  maneuvers  with 
vibration  suppression  for  a  flexible  space  vehicle  (Junkins, 
1989,  pp.  1-4).  The  central  hub  pivots  on  a  ceramic  bearing 
in  the  horizontal  plane,  and  four  identical  cantile'-’ered 
flexible  appendages  with  endpoint  masses  protrude  from  this 
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base.  A  reaction  wheel  is  driven  by  a  DC  brushless  motor  to 
provide  a  system  torque.  Sensors  include  an  angle  encoder, 
strain  gauges  and  a  motor  tachometer.  Control  laws  were 
successfully  written  to  control  the  system;  however, 
nonlinear  f riction/stiction  occurred  as  a  result  of  the 
ceramic  bearing,  and  nonlinear  bending  moments  were  observed 
as  a  result  of  the  flexibility  which  hampered  accurate 
modeling  of  the  system  response. 

Old  Dominion  University  has  conducted  an  experiment  to 
investigate  the  slewing  of  flexible  structures  while 
simultaneously  suppressing  vibrational  motion  during  the 
maneuver  (Yang,  1989,  pp.l).  The  experimental  setup 
consists  of  a  trolley  on  which  is  mounted  a  long,  flexible 
beam.  The  trolley  is  driven  through  a  driver  pulley  and  a 
cable  transmission  system.  The  flexible  beam  is  rotated  in 
a  horizontal  plane  by  the  beam  motor.  Sensors  include 
strain  gages,  angular  potentiometers  and  a  tachometer  for 
the  beam  motor.  Again,  nonlinear  effects  were  discovered  in 
large  bending  deflections,  and  friction  of  the  cables  and 
trolley,  but  again  they  posed  no  problem  in  controlling  the 
system . 

Stanford  University  has  set  up  an  experiment  most  similar 
to  the  one  designed  here,  and  mutual  cooperation 
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between  Stanford  University  and  Naval  Postgraduate  School 
has  proved  invaluable  in  the  early  stages  of  design. 

A  two  link  manipulator  is  used  by  Stanford  University 
consisting  of  two  flexible  beams  pinned  at  the  shoulder,  and 
joined  at  the  elbow  by  a  limited  angle  torquer  motor 
(Oakley , 1988 ,  pp.  1-4).  Rotary  variable  differential 
transformers  measure  joint  angles,  and  a  CCD  television 
camera  tracks  endpoint  position.  The  control  equations  were 
derived  using  an  assumed  modes  method,  and  the  experimental 
results  concurred  with  the  simulation.  Again,  despite 
nonlinear  effects,  accurate  modelling  was  achievable. 

At  the  Naval  Postgraduate  School,  analytical/ 
experimental  studies  have  been  started  this  year  on  the 
attitude  control  of  flexible  spacecraft. 

C.  EXPERIMENTAL  DESIGN 

The  design  of  the  experimental  set-up  requires  a 
compromise  between  simulation  realism  and  the  practicalities 
of  the  hardware.  A  design  is  needed  that  will  be  an 
accurate  model  for  understanding  spacecraft  control 
problems.  For  this  reason,  the  experiment  was  constrained 
to  rotate,  not  translate,  about  its  principal  axis.  Pitch 
motion  is  the  only  motion  being  investigated  and  controlled. 
In  space,  a  momentum  wheel  can  be  used  to  control  pitch  axis 
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motion  since  the  pitch  axis  is  the  only  axis  that  remains 
inertially  fixed  for  Earth-oriented  satellites. 

There  are  seven  phases  of  development  in  the  experimentc 1 
design.  This  phased  approach  to  building  the  experimental 
laboratory  allows  interim  goals  to  be  identified  and 
realized , 

1.  Phase  I 

In  the  initial  phase,  the  design  will  consist  of  a 
mainbody,  an  L-beam  flexible  arm,  and  a  reflector  (Figure 
1) ,  The  pitch  motion  wil^  be  controlled  by  a  motor  driven 
momentum  wheel.  There  will  be  an  angle  and  angul' r  rate 
sensor  on  the  mainbody.  The  flexible  arm  end  point  position 
will  be  measured.  Beam  mode-shapes  will  be  determined  bv 
strain  gages  and  accelerometers. 

The  mainbody  will  be  controlled  with  information  from 
the  ma ' nbody  only.  The  arm  sensors  will  be  utilized  for 
performance  measurements  and  mode  shape  analysis.  The 
reflector  will  not  be  controlled  during  this  phase. 

2.  Phase  II 

In  phase  II,  tne  central  body  will  be  controlled  with 
information  from  the  mainbody  rate  and  angular  position 
sensors  and  the  endpoint  position  sensor. 

3.  Phase  III 

A  limited-angle  stepper  motor  will  be  added  in  phase 
III  at  the  reflector-arm  joint  to  control  the  orientation  of 
the  reflector  with  respect  to  a  fixed  point  on  the  mainbody. 
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The  angular  position  of  the  reflector  with  respect  to  the 
end  point  will  be  measured.  The  central  body  will  be 
controlled,  as  above.  The  reflector  will  also  be  controlled 
in  both  open  loop,  closed  loop  with  feedback  from  sensors  on 
the  reflector,  and  closed  loop  with  feedback  from  all 
sensors . 

4.  Phase  IV 

Thrusters  will  be  added  to  the  mainbody  to  desaturate 
the  momentum  wheel. 

5.  Phase  V 

The  strain  gages  and  accelerometers  will  be  fed  into 
the  control  loop  for  controlling  the  reflector. 

6.  Phase  VI 

Liquid  tanks  will  be  added  to  the  system. 

7.  Phase  VII 

A  limited  angle  stepper  motor  and  angle  position 
sensor  will  be  added  at  the  L-joint  and  active  control  of 
the  flexible  arm  will  be  achieved.  Eventually,  robotics 
experiments  can  be  explored  by  fixing  the  mainbody. 
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D.  OBJECTIVES 


This  thesis  is  concerned  with  the  initial  design  of  the 
experiment  through  Phase  T.  There  are  three  main 
obj  ectives . 

1.  Experimental  Set-up 

The  preliminary  design  of  the  experimental  set-up 
includes  the  identification  of  system  requirements,  the 
resulting  specifications,  and  the  selection  of  sensors, 
actuators,  and  computer  systems  to  make  the  experiment 
operational  through  Phase  I. 

2.  Equations  of  Motion 

The  equations  of  motion  for  the  experimental  system 
must  be  written  to  determine  the  mass  and  stiffness 
matrices . 

3.  System  Identification 

The  system  characteristics  (the  natural  frequencies 
and  mode  shapes) ,  should  be  determined  by  computer  aided 
finite  eleinen*  :^nalysis  and  verified  by  experimental  modal 
analysis . 
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II.  EXPERIMENTAL  SETUP 


The  initial  design  of  the  experimental  set-up  consisted 
of  determining  the  required  specifications,  researching  the 
available  equipment,  and  selecting  the  appropriate  product. 
Whenever  possible,  the  set-up  was  designed  as  a  scaled-down 
version  of  existing  space-platforms.  The  specifications 
were  also  often  based  upon  real-world  considerations,  such 
as  acceptable  sizes  and  weights  for  a  laboratory  set-up. 

A.  GRANITE  TABLE 

Simulation  of  a  space  environment  requires  that  gravity 
and  friction  be  negligible.  A  large  granite  table  serves  as 
the  foundation  for  the  laboratory  set-up  because  precision 
flatness  and  smoothness  can  be  achieved.  The  mainbody  and 
flexible  arm  float  on  air  pads  that  skim  over  the  surface  of 
the  granite  table. 

The  ta^^le  is  6 '  X  8 '  X  10.5"  and  rests  on  a  castor  stand. 
There  are  set-screws  on  the  metal  stand  to  achieve  a  level 
platform.  The  surface  is  a  Laboratory  Grade  A  .001"  finish. 
The  stand  has  wheels  to  allow  mobility.  The  table  weighs 
7550  lbs. 
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B.  MAINBODY 


The  overall  size  of  the  mainbody  is  constrained  by  the 
dimensions  of  the  granite  table.  In  order  to  have  adequate 
room  to  maneuver  the  system  when  slewing,  a  30  inch  diameter 
was  chosen. 

The  base  of  the  mainbody  is  a  simple  7/8"  aluminum  disk 
with  a  15  inch  radius  (Figure  2) .  The  mainbody  was  machined 
at  Naval  Postgraduate  School  by  the  Aeronautical  and 
Astronautical  Engineering  department.  The  top  side  of  the 
mainbody  has  seventy  two  1/4  inch  #20  drilled  holes 
positioned  in  a  radial  pattern  to  allow  for  a  wide 
assortment  of  attachments  as  the  experiment  grows. 


10 


To  achieve  a  realistic  simulation  of  actual  spacecraft 
conditions,  the  moment  of  inertia  of  the  mainbody  is  assumed 
to  be  roughly  three  times  the  moment  of  inertia  of  the 
flexible  arm.  From  this  specification,  the  moment  of 
inertia  of  the  mainbody  can  be  calculated: 


^mainbody  ^  ^ arm  assembly'^ 

^mainbody  =  (3)  X  25,000  Ibw-inches^  (1) 

^mair^ody  =  75,000  IbW-jncheS^ 


The  total  mass  for  the  mainbody  can  be  calculated  from 
its  moment  of  inertia: 


T"  “  ^  "K  ^ 

■^mainbody  ~  ~^^mainbody^ inainbody 

^mainbody  =  ^67  IblTl 


Where  r„,.  =  15  inches 

ma  1 nbody 


(2) 


For  practicality  in  handling,  the  mass  of  the  entire 
mainbody  should  be  under  500  lb.  This  means  that  the 
inertia  ratio  of  flexible  assembly/mainbody  will  be  more  on 
the  order  of  2  to  1,  or  an  of  approximately  50,000  lb- 
in^.  This  gives  a  total  moment  of  inertia  for  the  entire 
system  of  approximately  75,000  Ib-in^. 

The  momentum  wheel  and  motor  assembly  will  add 
approximately  30  lb  to  the  system.  From  the  above 
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calculation,  it  is  obvious  that  additional  weights  must  be 
added  to  the  base  to  achieve  the  desired  moment  of  inertia. 

The  mass  of  the  base  of  the  mainbody  can  be  calculated  as 
follows : 


M, 


base 


base 


P  ^  ^ base^base 
=  59 . 376  Ihm 


(3) 


where p  =  .096  Ib-in^ 
^  base  =  inches 

h  base  =  inches 


From  the  above,  the  moment  of  inertia  of  the  base  of  the 
mainbody  can  be  calculated  by: 

1 

^ mb  ~  ^mb^mb  (4) 

^mb  ~  6,67  9.80  Ibm- inches'^ 


C.  FLEXIBLE  ARM 

The  flexible  arm  (Figure  3)  consists  of  two  assemblies  of 
aluminum  and  steel  construction  bolted  together  rigidly  in  a 
right  angle  elbow.  Mass  intensifiers  are  connected  to  the 
thin  (.16  inch)  aluminum  bar  to  increase  the  moment  of 
inertia  of  the  arm  assembly  without  significantly  increasing 
the  stiffness.  Each  mass  intensifier  (Figure  4)  consists  of 
two  parts  which  are  bolted  together  on  each  side  of  the 
aluminum.  The  tapered  design  of  the  mass  intensifier 
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reduces  the  contact  between  the  steel  and  the  aluminum. 
Minimal  contact  will  ensure  greater  flexibility. 


Figure  3  Flexible  Arm  Assembly. 
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Figure  4  Mass  Intensifiers . 


The  elbow  joint  consists  of  two  aluminum  right  angle 
braces  connected  to  a  circular  aluminum  disk.  The  end  joint 
has  an  aluminum  wedge  bolted  to  a  circular  aluminum  disk. 
This  wedge  shape  allows  free  motion  of  the  endpoint.  The 
end  of  the  flexible  arm  will  house  a  reflector  and  will  be 
the  location  for  end  point  sensing. 

The  total  mass  of  the  flexible  arm  assembly  can  be 
calculated  as  follows: 

mass  of  aluminum  rod^  =  .411  Ibm 
mass  of  aluminum  rodg  =  .380  Ibm 
mass  of  each  intensifier  pair  =  1.056  Ibm 

Total  mass  of  arm^^  =  9.66  Ibm 
Total  mass  of  arm^  =  8.59  Ibm 

mass  of  elbow  joint  =  .244  Ibm 
mass  of  end  joint  =  .183  Ibm 

Total  mass  of  flexible  assembly  =  18.68  Ibm 
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The  center  of  mass  of  the  flexible  arm  assembly  is: 


Centex  of  mass  of  arm. 


CM. 


12 . 96  inches  1 


Center  of  mass  of  armg  = 


E 


Af. 


CM. 


arnin 


=  11.66  Inches  j 


Center  of  mass  of  arm  assembly  ' 

assembly 
^y)  assembly  ~ 


cm 

M 


CM^r^  (x)  +  M^^^  CM^rrr,  ( ^) 

ariBj.  aztDx  Ar/n*  ' 


Maza:.CMaza,.  (v)  *Kza.CM,,„^  (y> 


^a  rm^  "^^a  zrrtg 


=  19.29  inches 
CM f.y)  assembly  =  5.55  inches 


(6) 


The  moment  of  inertia  of  each  arm  can  be  calculated  about 
the  center  of  mass  of  each  arm.  Each  mass  intensifier  and 
the  elbow  and  end  joint  braces  are  computed  as  point  masses. 
The  moment  of  inertia  for  the  flexible  arm  assembly  can  be 
calculated  about  the  origin  O  by  the  parallel  axis  theorem. 
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)  =  416.10  Ibm-inches^ 

dcM  ^  -  324.03  Ibirt-inches" 


M. 


where  d  -  (15+12.96)  =  27.96  inches 


(J,^_  )  =  7,967.9  2  Ibm-inches^ 


(7) 


( ) 


^  ^armg 


where  d  =  v^(15  +  27 . 96  )  ^+ ( 11 . 66  )  ^  =  42.85  inches 


(J...  )  =  16,096.32  Ibm-inches^ 
*0 


■^arm  assemhiy^  ^ 

^a^assemli.ve  =  24,064.24  Ibm-inches^ 


D.  MOMENTUM  WHEEL 

A  motor  driven  momentum  wheel  will  be  used  to  apply  a 
torque  (the  change  in  the  angular  momentum)  for  slewing  the 
ma inbody . 

For  preliminary  calculations,  a  one  foot  diameter 
aluminum  disk  with  one  inch  thickness  can  be  used  for  the 
dimensions  of  the  momentum  wheel. 
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The  mass  of  the  momentum  wheel  can  be  calculated  as: 


P  ^  ^mw  ^mw 

Where  r^^ 

=6.00  inches 

^mw 

=1.0  inches 

P 

.  .096 

inches 

=  10.857  Ibm 

(8) 


The  moment  of  inertia  for  the  momentum  wheel  is: 


I„ 


where 


I. 

I. 


6.00  inches 
10.857  Ibw 
195.43  Ibw-inches^ 
8.094  oz-inches-sec^ 


(9) 


E.  TORQUE  MOTOR 

The  momentum  wheel  must  be  driven  by  a  motor  that  will 
deliver  the  required  change  in  angular  momentum  to  the 
system.  A  specification  for  the  system  design,  based  upon 
existing  space  platforms,  is  that  the  mainbody  accelerate  at 
the  rate  of  12  degrees  per  second^. 
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From  this  specification,  the  torque  of  the  motor  can  be 


determined  to  be: 


^syscem^syscen 


(3000oz-irches-sec^)  (.209 


rad 

sec^ 

a 


]  =  (8 . 094oz-inches-sec^)  (a  ) 


mw 


=  80 


rad 

sec^ 


(10) 


Torque  of  momentum  wheel  =  a^. 

=  6  50  oz-inch 


From  the  above  calculation,  the  momentum  wheel  must  apply 
a  650  oz-inch  (4.59  N-m)  torque  to  provide  the  desired  slew 
rate  for  the  mainbody. 

The  motor  selected  was  a  DC  servo-disc  motor.  In  the 
servo-disc  motor,  the  ironless  armature  is  constructed  from 
several  layers  of  copper  conductors  in  a  flat-disc 
configuration.  This  allows  for  low  inertia  and  fast 
acceleration.  Additionally,  the  smaller  size  can  be  easily 
mounted  on  the  mainbody.  Table  I  outlines  the  motor 
characteristic . 


18 


TABLE  I  PMI  SERVO-DISC  MOTOR 

Motor  Performance 

Units 

TR16M4C 

Peak  Torque 

oz-in 

5307.2 

Continuous  Stall  Torque 

oz-in 

498 . 4 

Peak  Current 

Amps 

100.8 

Peak  Acceleration  w/o  Load 

KRADS/sec^ 

63 . 2 

Horsepower 

HP 

1.4 

Torque 

oz-in 

Speed 

RPM 

3000 

Power  Output 

Watts 

1049 . 5 

Terminal  Voltage 

Volts 

128 . 7 

Current 

Amps 

9.55 

Torque  Constant  (KT) 

oz-in/ 

52.77 

Motor  Weight 

Lb 

17.50 

Moment  of  Inertia 

oz-in-s^ 

.  084 

Cost  (includes  tachometer) 

ea 

$1229 . 0 

The  motor  will  be  aligned  in  a  vertical  position  on  top 
of  the  momentum  wheel  to  deliver  a  torque  about  the  vertical 
(pitch)  axis.  The  motor  housing  includes  an  attachment 
shaft  for  fitting  with  the  momentum  wheel. 
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F.  SENSORS 


Initially,  the  angular  position  and  the  angular  rate  of 
the  mainbody  and  the  speed  of  the  momentum  wheel  will  be 
sensed.  At  later  time,  the  angular  position  of  the 
reflector  with  respect  to  the  endpoint,  and  the  position  of 
the  endpoint  with  respect  to  the  mainbody,  will  be  measui  d. 

The  angular  rate  of  the  torquer  mooor  and  momentum  wheel 
will  be  sensed  by  a  tachometer  which  is  attached  to  the 
motor.  Table  II  outlines  the  tachometer  characteristics. 


TABLE  II  PMI  TACHOMETER 

Tachometer  Performance 

Units 

JR16M4CH 

Maximum  Speed 

RPM 

4000 

Bidirectional  Tolerance 

%V 

±1  •  5 

Tachometer  Weight 

lb 

2 . 2 

Moment  of  Inertia 

oz- 

.  005 

in-s^ 

The  angular  rate  of  the  mainbody  will  be  monitored  by  a 
rate  sensing  gyroscope.  The  gyroscope  produces  an  analog 
output  voltage  proportional  to  the  angular  velocity  about 
the  sensing  axis  ''in  this  case,  the  Z  axis  at  the  origin). 
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When  a  rotation  exists,  the  Coriolis  forces  transfer 
momentum  to  the  perpendicular  plane  and  cause  a  bending  of 
the  solid  state  sensing  elements.  These  elements  will 
vibrate  180  degrees  out  of  phase  with  angular  motion  but  in 
phase  when  linear  acceleration  or  vibration  occurs.  The 
integrated  electronics  produce  an  angular  motion  response 
output  with  no  linear  motion  incorporated.  A  single  axis 
gyroscope  with  a  range  of  +30  degrees  per  second  was 
selected.  Table  III  shows  the  rate  gyroscope 
characteristics . 


TABLE  III  WATSON 

RATE  GYROSCOPE 

Rate  Gyroscope 

Units 

ARS-C121-1A 

Output 

VDC 

0  at  zero  angular  rate 
+10  at  full  scale 
angular  rate 

Power  Supply 

VDC 

+  15  +  5%  20  mA 
maximum 

Range 

degree 

/sec 

+  30 

Weight 

Oz 

4 

Cost 

ea 

$  802.00 

Angular  position  for  both  the  mainbody  and  the  endpoint 
reflector  will  be  measured  by  rotary  variable  differential 
transformers  (RVDTs) .  RVDTs  measure  the  angular 
displacement  of  rotating  elements  by  producing  a  voltage 
whose  magnitude  varies  linearly  with  the  angular  position  of 
the  shaft.  The  RVDT  is  constructed  on  precision  ball 


bearings  to  minimize  friction  and  noise.  The  best  operating 
range  of  a  RVDT  is  +  40  degrees,  which  is  well  within  the 
confines  of  this  experiment.  The  characteristics  of  the 
RVDTs  purchased  are  in  Table  IV. 


TABLE  IV  SCHAEVITZ  ROTARY 
VARIABLE  DIFFERENTIAL  TRANSFORMER 


RVDTs 

Units 

R30D 

Linearity 

%  of 

+  3  0°  :  .25 

range 

+  40°  :  .50 

+  60°  :  1.5 

Moment  of  Inertia 

Oz-in- 

8.5  X  10’'^ 

sec^ 

Maximum  Torque 

Inch-oz 

.  019 

Maximum  Load 

Lb 

Radial  =  8 

Axial  =  10 

Weight 

Oz 

1.87 

Cost 

ea 

$438 . 00 

Endpoint  sensing  has  yet  to  be  determined.  The 
accuracies  required  (+  .1  mm)  over  the  large  span  of  control 
(+  5  cm)  make  most  conventional  endpoint  sensing  techniques 
inadequate.  CCD  cameras  would  provide  the  required  coverage 
area,  however  they  would  give  only  1  mm  accuracies  and  the 


signal  processing  equipment  required  to  integrate  CCD 
cameras  into  the  control  loop  would  be  expansive.  Greater 
accuracies  can  easily  be  achieved  by  off-the-shelf 
capacitance  sensors,  however  the  coverage  area  would  be  much 
smaller  than  that  required  by  this  experiment. 

A  promising  devise  for  lateral  position  sensing  is  the 
electro-optical  photodetector.  Using  a  fixed-position  laser 
or  LED  light  source,  a  mirror  system  and  a  photodiode 
receiver  collect  the  reflected  light  and  provide  an  output 
signal  proportional  to  the  position. 

Standa^-d  off-the-shelf  units  can  be  purchased  through 
several  sources.  The  unit  would  consist  of  an  analog 
optical  position  monitor  (approximate  cost  $3,500)  and  a 
detector/cable  assembly  (approximate  cost  $400)  . 

G.  AIR  PADS 

The  mainbody  and  the  flexible  arm  will  be  supported  by 
five  aluminum  air  pads,  each  capable  of  supporting  a  60  lb 
load.  Each  pad  is  fed  by  80  PSI  of  air,  which  is  regulated 
by  two  regulators.  One  regulator  feeds  the  three  mainbody 
air  pads,  and  one  regulator  feeds  the  elbow  and  endpoint  air 
pads . 

On  the  mainbody,  three  1.5  inch  diameter  mounting  holes 
are  located  120  degrees  apart  three  inches  from  the  outer 
edge.  The  air  feed  is  on  the  side  of  each  pad.  A  top  brace 
will  be  fashioned  and  mounted  to  the  mainbody  to  hold  the 
air  pads  in  place.  Three  air  pads  for  the  mainbody  have 
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been  chosen  as  the  optimum  number  for  controlling  the 
mass.  An  odd  number  of  pads  decreases  the  likelihood  of 
leverage  effects  due  to  uneven  load  distributions.  A  150 
psi  regulator  will  control  all  three  mainbody  air  pads.  The 
150  psi  regulator  was  chosen  to  allow  upgrades  to  larger  air 
pads  at  a  later  date. 

H.  AIR  BEARING 

The  entire  system  will  be  constrained  to  rotational 
motion  only  by  an  air  bearing  mounted  to  the  center  of  the 
mainbrtdy  from  above.  Translational  motion  will  be 
constrained  in  the  initial  phase  of  the  experiment  so  that 
pitch  effects  can  be  isolated.  At  a  later  time,  the  air 
bearing  can  be  removed  so  that  robotics  experiments  can  be 
performed.  The  air  bearing  (Figure  5)  consists  of  an  H-type 
bearing  which  is  capable  of  carrying  both  a  radial  and  an 
axial  load.  As  the  mass  of  the  mainbody  is  increased,  the 
air  bearing  can  be  used  to  support  some  of  the  load.  A 
mounting  bracket  will  be  attached  to  the  side  of  the  granite 
table  and  will  be  connected  to  the  mainbody  from  above. 
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Figure  5  Air  Bearing. 


I.  COMPUTER  SYSTEM 

The  AC-100  computer  system,  manufactured  by  Integrated 
Systems,  Inc.,  has  been  purchased  to  provide  data 
acquisition  and  results  analysis  for  the  experimental  set¬ 
up.  The  AC-100  is  an  automatically-programmable  real-time 
control  system  that  will  allow  graphical  specification  and 
simulation  of  real-time  systems.  This  will  allow  the 
experimental  set-up  to  be  fully  monitored  real-time  from  the 
workstation. 

The  AC-100  system  consists  of  three  parts;  the 
workstation,  the  AC-100  software,  and  the  AC-100  controller. 
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The  workstation  includes  a  DEC  3100  computer,  VMS  operating 
system,  8  MB  RAM  memory,  two  105  MB  hard  disks,  a  19  inch 
color  monitor  and  an  Ethernet  interface.  PWS003  modeling, 
simulation  and  implementation  software  are  combined  with  a 
VAX  C  and  VAX  Fortran  compiler.  The  controller  has  16 
channel  inputs  and  10  channel  outputs  in  addition  to  32 
parallel  digital  inputs/outputs,  eight  encoder  inputs,  and 
one  event  trigger  input. 
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III.  EXPERIMENTAL  ANALYSIS 


A.  MODAL  ANALYSIS 

Experimental  modal  testing  can  be  conducted  to  determine 
the  nature  of  the  vibration  response  and  to  verify  the 
analytical  models.  The  Hewlett-Packard  VISTA  Data 
Acquisition  Package  was  used  to  collect  information 
concerning  the  time  and  frequency  response  of  the  flexible 
arm  assembly  in  a  free-free  mode. 

The  frequency  response  function  method  of  modal  testing 
requires  that  the  input  excitation  and  the  output  response 
be  measured  simultanously  to  determine  the  system  function 
(Figure  6) . 


The  VISTA  program  utilizes  a  dynamic  signal  analyzer, 
which  is  a  Fourier  transform-based  instrument,  to  process 
the  measurements  of  the  structures  frequency  response 
(Hewlett-  Packard,  1986,  pp.  11). 
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The  input  used  was  a  single-point  excitation  impact 
hammer.  The  hammer  is  the  most  simple  way  to  excite  a 
structure  into  vibration.  The  exciter  consists  of  an 
impactor,  with  several  different  tips  and  heads  which  will 
allow  some  flexibility  in  setting  the  force  level  ranges  for 
testing  different  objects  (Ewins,  1984,  pp.  102-104).  The 
hammer  contains  an  ICP  quartz  force  sensor  mounted  on  the 
striking  end  of  the  hammer  (PCB  Piezotronics,  1983,  pp.  1). 
The  impact  force  is  transferred  via  the  sensing  element  into 
an  electrical  signal  which  can  be  evaluated. 

The  hammer  impulse  consists  of  a  nearly  constant  force 
over  a  wide  frequency  range  and  therefore  can  excite  all 
resonances  in  that  range.  The  size  of  the  hammer,  along 
with  the  hammer  tip  material  and  velocity,  determine  the 
amplitude  of  the  force  impulse.  The  frequency  content  of 
the  energy  applied  to  the  arm  is  a  function  of  the  stiffness 
of  the  contacting  surfaces  and  the  mass  of  the  hammer.  The 
hammer  tip  selected  affects  the  force  impulse,  and  therefore 
the  frequency  response.  A  hard,  steel  tip  was  selected  to 
deliver  a  short  pulse  duration  and  higher  frequency  content. 

The  flexible  arm  assembly  was  supported  from  above  by 
small  soft  elastic  cords  in  an  attempt  to  simulate  the  free- 
free  mode  (Figure  7) .  An  ICP  accelerometer  was  positioned 
in  the  XY  plane  at  the  end  of  each  arm.  When  the 
accelerometer  vibrates,  an  internal  mass  applies  a  force  to 
the  crystal  element  which  is  proportional  to  the 
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acceleration.  Using  Newton's  Law,  the  resulting  force 
output  can  be  measured. 


Figure  7  Flexible  Arm  Assembly  Support. 


The  assembly  was  struck  with  the  hammer  in  various  places 
throughout  the  assembly.  The  resulting  frequency  response 
showed  consistent  modes  of  varying  amplitudes,  depending 
upon  where  the  structure  was  struck. 

Figure  8  shows  the  time  and  frequency  spectrum  for  the 
impact  hammer  as  struck  about  midway  on  arm  A.  Figure  9  and 
Figure  10  show  the  time  and  frequency  measurements  for  each 
arm  as  a  result  of  the  impact. 
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T  1  rna 

R;  Time  Record 


Figure  8  impact  Hammer  Time  and  Frequency  Measurements. 
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Figure  9  Arm  A  Time  and  Frequency  Response. 


Figure  10  Arm  B  Time  and  Frequency  Measurement. 


31 


The  finite  element  analysis  was  evaluated  for  the  free- 


free  condition,  and  the  mode  frequencies  were  computed. 

When  compared  with  the  experimental  results,  a  fairly  good 
correlation  was  found  between  the  two,  with  the  exception  of 
the  first  mode.  Table  V  shows  the  difference  between  the 
frequencies  obtained  from  the  GIFTS  program  and  the 
experimental  values. 


TABLE  V  FINITE  ELEMENT  ANALYSIS 

VS.  EXPERIMENTAL  RESULTS 

MODE 

GIFTS 

Hz 

EXP.AVG 

Hz 

DIFFER¬ 

ENCE 

Hz 

1 

.  548 

1.125 

-.577 

2 

2.087 

2.249 

-.162 

3 

3.057 

3.375 

-.318 

4 

6.483 

6.093 

.  390 

5 

8.133 

8.312 

-.179 

6 

13 .424 

12.313 

1.111 

7 

15.667 

15.125 

.  542 
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IV.  THEORETICAL  ANALYSIS 


A.  EQUATIONS  OF  MOTION 

The  system  equations  of  motion  can  be  derived  as  a 
mathematical  model  for  the  exoerimental  set-up.  This  model 
is  based  upon  given  physical  laws  which  govern  the  system 
behavior,  such  as  force-deformation  relationships  and  the 
laws  of  motion.  The  equations  of  motion  wil?  identify  the 
mass  and  stiffness  system  parameters. 

In  a  "lumped-parameter”  system,  the  equations  are  a 
function  of  time  alone,  which  are  ordinary  differential 
equations  and  are  easy  to  solve  (Meirovitch,  1990,  pp.  45- 
53).  In  the  "distributed-system-parau.eter"  system,  the 
equations  are  functions  of  both  spatial  coordinates  and 
time,  and  are  therefore  partial  differential  equatic  .s  which 
can  not  be  easily  solved.  Therefore,  a  discret  model  must 
be  devised  to  model  the  distributed  system. 

A  flexible  system  is  a  non-linear  system.  In  other 
words,  the  response  of  the  experimental  set-up  to  different 
excitations  can  not  be  obtained  separately  and  combined 
linearly.  Nonlinearity  may  present  complications  in 
modelling  and  control.  A  careful  ccmparison  of  the 
analytical  verses  experimental  results  will  identify 
nonlinear  variations. 
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The  equations  of  motion  for  a  flexible  structure  can  be 
written  as  a  function  of  spatial  variables,  which  can  be 
assumed  cjs  known,  and  unknown  time  variables.  This  will 
result  in  a  discrete  model  containing  ordinary  differential 
equations.  The  solution  can  be  represented  as  spatial 
eigenfunctions  multiplied  by  time-dependant  generalized 
coordinates . 

Lagrange's  equation  facilitates  writing  the  equations  of 
motion  for  a  flexible  multiple  degree  of  freedom  system. 
Additionally,  Lagrange's  technique  produces  an  entirely 
scalar  process  using  the  quantities  of  potential  energy, 
kinetic  energy,  and  work  in  terms  of  a  generalized 
coordinate  system  (Meirovitch,  1967,  pp. 47-59). 

The  Lagrange  Operator  is: 


d  ,  dL 
d  C  dq^ ' 


(It) 


where  L  is  the  Lagrangian:  T  -  V 
T  =  System  kinetic  energy 
V  =  System  potential  energy 
q-  =  System  generalized  coordinates 
(/>  =  SyF.:e  .  generalized  force 
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By  using  this  method  of  analysis,  the  complicated  dynamic 
problem  can  be  reduced  to  a  scalar  integral  that  does  not 
depend  on  the  coordinates  used. 

The  kinetic  energy  of  a  distributed  system  is: 

T=f^p(V^)^  ds  (12) 

s 


where  p  is  the  material  density,  is  the  velocity  of  a 
generic  point  of  the  system  in  inertial  frame 


The  origin  O  for  the  experimental  set-up  was  chosen  to  be 
the  center  of  the  mainbody.  This  is  the  point  where  the 
mainbody  is  attached  to  the  air  bearing  (Figure  11) .  This 
origin  O  is  therefore  fixed  in  the  reference  frame  of  the 
table . 

The  value  of  can  be  determined  from  the  following 
formula : 


where  o)  =  D  ^ 


(13) 


35 


36 


The  position  vector  from  the  origin  O  to  any  point  on  the 
arm  can  be  described  by: 

?o-a  =  X  i  +  j 
^o-b  =  (R+Li+W,)  i  +  +  y)  j 

where  =  length  of  Arm  A 
R  =  radius  of  the  mainbody 
Wgg  =  the  vertical  endpoint  displacement  of  Arm  A 

thus 

w  X  fo.3  =  X  [  ixi+W^j)  ]  = 

o3  X  =  6ic  X  [  (R+Li  +  f^)  1  +  (y+Wag)  j]  (14) 

=  1  +  (i?+Li+f4)6  j 

Using  the  function  Wg(x,y,t)  to  represent  the  position 
of  any  arbitrary  point  on  beam  A,  and  to  represent 

the  position  of  any  arbitrary  point  on  beam  B: 

W,  =tira,(x)g,(t)+ilr^2(x)g2(t) 

+tb2(y)g2(«^)  mc, 

=  +1^32  ( L2 )  g^  (  t) 

=Pi  g^  ( t)  +P2  Qz  ( t) 

where  t  is  a  function  of  spatial  variables  only  and 
q(t)  is  the  time  dependant  general  coordinate 
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The  subscripts  1  and  2  refer  to  the  first  and  second 
inodes  of  the  structure.  Only  the  first  two  inodes  have  been 
incorporated  in  the  spatial  equations,  but  others  could 
easily  be  added  at  a  later  time.  The  displacement  in  X 
direction  for  Arm  A  has  been  ignored  as  insignificant 
compared  to  the  other  displacements. 
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Substituting  this  new  variable  into  the  equation,  the 


kinetic  energy  of  the  system  can  be  calculated  to  be: 


KEsystein=  —  I^ 


R*L, 


^0  R  ^0 

R*Li 

[pJpiPjdy+p  I  4r^_j4ra2<*<^-*-p/4'ii^i,2C?y] 

OR  0 

L2  R*L^  L2 

[-|/p2C?y+^  /  'ifl2dx+-^j}\il2dy] 

0  R  ^  0 

R*Li  L2 

R  ^  0 

R*L^  L2 

[PP1P2L2  +  P  I  taj4'a2^+p/’ifijllffc2C?y] 

R  0 

L2 

[p  (i?+Li)L2pi+P  /  X}\faidx-pfYy^j,idy] 

R  0 

R*Li  L2 

[pp2L2+-|  J  tLdx+^JtLciy] 

i?  ^0 

i?  * 

[p  (i?i+L^)L2P2+p  f  x^^zdx-pfy^bzdy'i 

R  0 

R*Li  Lj 

[-^  j  x^dx+-^  Jy^+ (i?+Li)  ^dy] 

[pj’(i?+Lj  i|;^^dy]  q^B^+[pj  (i?+Li) 

0  0 

J^2  ^ 

tp/ p2'l'fci^y]  diQ2^  -  [p[ Pztfc^dy] 

0  0 

[p|p2rdy]  g2^^  [p/Pii"c?y] 

0  0 
L2  1*2 

ip/p,*  t2<iy]  gi^zP  -  tp/PitfciCfyl 


e^+ 

gfe^- 

<7iQ-2p"  + 

<7iP  + 

g2p  + 

g20^  - 
gzdjP  + 

g^e^  - 

<7igi0 


(16) 
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The  potential  energy  of  the  system  is  the  flexural  strain 
energy.  To  calculate  the  potential  energy  of  each  beam,  the 
following  equation  can  be  used: 


y-l  /  EI( 


dx^ 


)  ^dx  + 


1 

2 


dX 

dy^ 


)  ^dy 


(12) 


Substituting  the  expressions  for  W^,  and  into  this 

equation,  the  total  potential  energy  of  the  system  is 
calculated. 

The  rotary  inertia  effect  and  the  axial  deformation 
effect  on  the  potential  energy  have  been  ignored.  The  cross- 
sectional  dimensions  are  small  compared  to  the  length  of  the 
bar. 
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The  Lagrangian  can  now  be  written  in  the  following  form 


(18) 


R  ^0 

R*L^  1/2 

[■|/P?c?y+-|  f  qf6^  + 

0  J?  0 

Z/^  R*L>y  Z/2 

[pJpjPjdy+p  f  i|raj»|»a2cbf+p|4»i,j’l»i,2<^y]  didjP"  ^ 

0  j;  0 

X.J  x,2 

[-|/P2c?y+-|  /  i|'Lci^+ /i|»Lc(y]  + 

0  ^  0 

^♦Lj  L2 

+  /  ’l'L‘i^+-|j‘i(r|idy]  + 

L2 

CPPlP2-t'2+P  /  'l'flJ»l»a2CiX+p/ti>2'l'i,2'^y]  + 

R  0 

R*L,  Lj 

[p(/?+Lj)L2Pi  +  p  J  Xilr^jdx-p Jyilr^jdy]  dj6  + 

J?  0 

R*I^  L3 

[-|P2i2  +  -|  f  ^l2dx*-^j-^l2dy]  + 

R  ^0 

Li 

[p  (i?  +  Ll)L2P2  +  P  I  X^^2‘^-pfy^b2^y'i  <r2^  * 

R  0 

[p|(i?+Li)i|rj^jdy]  gi6^+ [pj  (R+Lj  il/^^dy]  - 

0  0 

i-a  Lj 

tp/twP2cfy]  gzP  -  p/P2»l>i,2dy]  q.d[,6  + 

0  0 

[p/P2i'cfy]  <72®^  +  [p/Pi5^c?y]  diP^  - 

0  0 

Li  Lj 

tp/H'fciPi'^y]  gi<?2P  -  [p/’kwPiC?y]  gi^iP  - 
0  0 

R*L.  Li 


J  dx^  dx^  J 


)dy]  gjgj 


jj  —  0  <^y^ 
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Performing  the  Lagrangian  Operator,  the  resulting  matrix 
is  as  follows: 

Mil  M12  M13  Kll  K12  K13  Qi  0 

M21  M2 2  M2 3  Qz  +  ^21  K22  K23  ^2  +  ^2  =  0 

M31  M32  M33  0  K31  K32  K33  0 

where  : 


Mil 

Ml  2 

Ml  3 

M21 

M2  2 

M2  3 

M31 

M3  2 

M3  3 


(20) 
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Jj  dx^  J  dy2 

Ei?:flii^cix*Ei1^  ^ 

•jj  dx^  dx^  J  dy^  dy^ 

^dx.Ejf -£!iM  ^  dy 

•jj  dx^  dx^  J  dy2  dy^ 

£’J  r  ( )  ^dx+EI  f  (  — ^  )  2dy 

K23  =  K31  =  ii:32  =  K33  =  0 


Lj 

Vi  =  -lpfpldy->-p  I  ^iicbc-t-p^i^l^dy]  - 

0  R  0 

R+I..  Lj 

[p/PiPjCfy+p  /  - 

OR  0 

^2  J^2 

[pj  (i?+Li)  ilr^jdy- tp^Pi ydy]  6^  - 

0  0 

tp/p2’i»i.idy-p/piT|;j,^dy]  - 

0  0 

Ip f Pit i,i<^y^  ^26 
0 

tp/Pi4’w<^y]  <3^16 
0 

I/j  i?*Lj  1^ 

[p/PiP2dy+p  I  ’I'aitarCLx'-^p/twth^rfy]  diP"  - 

OR  0 

L2  R*Z^  Lj 

[pfpidy*p  J  llllzdx■^pJ1|^l2dy]  - 

0  ft  0 

L,  L, 

[pJP2ydy+pJ(i?+I,i)  i|r^^dy]  P^  - 

0  0 

t2pJP2i|fi,2dy]  (JjP  - 

0 

fp/Pi’l'wC?y"^p/'l'wp2c?y]  - 

0  0 

ip/Pi’i'wcfyi  <7iP  - 
0 

[p/P2t/>2Cfy^  ^2® 
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—i  —1  < 

V3  =  [2pj'pldy+2p  J  4r^idx+2pji|fwdy]  + 

0  A  0 

R*L2 

[2pj'pjdy+2p  f  }tiljdx*2pfililjdy]  ^1(^16  + 

0  J?  0 

Zpj  ^*Xh  i»2 

[2pj'PiP2dy+2p  J  i|r^i4'*2dx+2pji|i^_,\(r^^dy]  + 

0  0 

[2pJ(i?+I.i)i|j^jdy+2pj‘Piydy]  4ri6  + 

0  0 

[2pJ(i?+Lj  4ii2dy+p|P2ydy]  <^26  + 

0  0 

Z»2  ^♦X.j  Z»2 

[2p|piP2dy>2p  /  4'a2’lf«2dx-+2p|iltijT|ri,^dy]  gi(3r3e  4 

0  ^  0 

^  I, 

[p/Pirfy+p  /  tLci>f+p/4'Lc?y]  + 

0  /e  0 

R*Li-i  Lt2 

[2p|p,p2dy+2p  /  <l»a2tai<dy+2p/^liiji|»fc^dy]  g^gjP  f 
0  0  0 

X<2  ^♦Xn  X>2 

[pjpfdy+p  J  i|»Lcix+p J+lidy]  gffi  + 

0  «  > 

I,  X, 

[2p|(i?-^Li)i|>3,jdy+2pJP3ydy]  gjP  + 

0  0 

1-2  ^2  ^ 

[2p|(i?+I,i)\lfjb2C?y+2p|P2ydy]  g^e  -  [plPj’I'wCfy]  (Jigz  - 
0  0  0 

L,  L2 

[p/P2tbi<^y-p/Pitb2cfy]  (jr^gj  -[pJPzt^^dy]  - 

0  0  0 

[p/P2’l'w‘^y]  '72<3^2  +tp/Pl’l'b2‘*'5  <7l<?l  - 

0  0 

[p/Pi’i'w'^y]  -  [p/+t.jPiC?y]  <?i^i 
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The  spatial  variables  can  be  determined  in  several  ways. 
The  assumed-method  (Kane,  1983,  pp.  318-319)  makes  the 
approximation  that  the  deflection  of  each  flexible  link  can 
be  expressed  as  a  sum  of  a  finite  number  of  spatial  mode 
shapes  that  satisfy  the  given  boundary  conditions.  For  a 
uniform  cantilever  free-free  or  fixed-free  beam,  the 
components  of  the  spatial  variables  have  been  calculated  by 
numerous  sources. 

The  experimental  set-up,  however,  is  not  a  simple 
cantilever/f ree-f ree  continuous  beam  construction  and 
therefore  would  not  conform  to  the  above  simplification. 
Therefore,  a  finite  element  analysis  of  the  flexible  arm  was 
conducted  to  determine  the  first  two  mode  shapes.  The 
finite  element  analysis  program  GIFTS  was  used  to  create  a 
model  of  a  continuous  aluminum  bar.  Point  mass  loads  were 
introduced  at  the  position  of  the  mass  intensifiers  and  the 
aluminum  connecting  elbow  and  end  braces.  The  structure  was 
divided  into  forty  nodes. 

For  a  fixed-free  representation,  the  position 
was  supported  in  the  finite  element  analysis  model.  From 
the  GIFTS  analysis,  unit  deflections  for  each  mode  shape 
could  be  determined.  Only  the  first  two  modes  were 
analyzed.  Each  unit  was  arbitrarily  determined  to  be  one 
inch.  By  adding  the  unit  deflection  to  the  original 
position,  a  polynomial  could  be  fitted  for  each  mode  shape 
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T 


(Figui.es  12  &  13)  .  Therefore,  the  spatial  variables  are  as 
follows: 


=  -  . 0057 

- 

.0345 

X  + 

.  0679 

'I'm 

=  .0009 

y2 

.3233 

Y  + 

27.72 

'I'a2 

=  -.0025 

+ 

.2462 

X  - 

.4822 

'J'm 

=  .0074 

y2  ^ 

.1958 

Y  + 

25 . 17 

The  vertical  deflection  of  Arm  A  endpoint  can  be 
determined  by  evaluating  this  function  to  be: 


=  -4.77 
-  4.26 


RRI  I  B 

psini 

RRM  1 

MODE  ] 
RRM  B 

MODE  ] 
RRM  R 


Figure  12  MODE  1  FINITE  ELEMENT  ANALYSIS 

VS.  POLYNOMIAL  FIT. 


(24) 


Figure  13  MODE  2  FINITE  ELEMENT  ANALYSIS 

VS.  POLYNOMIAL  FIT. 


Performing  the  assigned  integration  results  in  the 
solution  of  the  matrix  variables  as  foliowu: 


Mil  =  202.84 
il2  =  -44.65 
M13  =-316.25 
M21  =  -44 .65 

M22  =  179.05  (25) 

M23  =  170.81 

M31  =-316.25 

M32  =  170.81 

M33  =  +  168.14 
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=  -  153.18  -18.57  - 

398.68  6^  +  39.84  q^^  - 
19.72  q^Q  -  22.08 

=  -155.7  8  q^6‘^  -  52.50  q^B^  - 

363.07  6^  -  35.94 
.4  4r^6  +  20.12  gi6  - 
17 .97  gjQ 

^3  =  305.00  gj^j^  +  306.36  (26) 

+  13.50  (Jr^gjO  +  13.50  + 

71.0.32  <^36  +711.11  ^26  + 

153.18  ql6  +  u.3.50  g3g26  + 

152 . 50  g^0  +  793.90  g^Q 
711.11  g20  -  19.72  g3g2  - 

39.84  qiQ2  -  17.97  - 

17.97  g^g,  -  22.08  g^g^  + 

22.08  g^  +  22  .  OSg^g, 


Kll  =  .71 

K12  =  .435 

K21  =  .435 

/C22  =  1.25 


(27) 
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B.  NATURAL  FREQUENCY 

When  placed  in  motion,  the  oscillation  of  a  system  that 
takes  place  will  be  at  the  natural  frequency  u^,  which  is  a 
property  of  the  system.  Damping  in  small  amounts  will  have 
little  effect  on  the  natural  frequency  and  has  been 
neglected  in  these  preliminary  calculations.  For  a  single 
degree  of  freedom  linear  system,  the  natural  frequency  can 
be  determined  from  the  equations  of  motion  by: 

=  I  (28) 

M 

where  K  is  the  stiffness  matrix 
M  is  the  mass  matrix 
For  the  experimental  set-up; 

E  =  10  X  10*  Ibm/in"^ 

I  =  2.0345  X  10'^  in^ 

=  .381  Ibm/in 

Ignoring  the  nonlinear  terms,  solving  the  equation  for 

=  .  0460  Hz 
(*>2  =  •  1004  Hz 
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V. 


CONCLUSIONS  AND  RECOMMENDATIONS 


Based  on  this  work,  the  following  conclusions  are 
reached . 

A.  SYSTEM  DESIGN 

The  experimental  set-up  for  Phase  I  study  of  spacecraft 
control/structural  interaction  has  been  designed.  All  the 
actuators  and  sensors  have  been  selected  except  the  end¬ 
point  displacement  sensing  of  the  arm.  The  mainbody  and  the 
flexible  arm  have  been  fabricated. 

B.  THEORETICAL  ANALYSIS 

The  equations  of  motion  for  the  experimental  model  have 
been  derived  and  natural  frequencies  determined. 

C.  EXPERIMENTAL  ANALYSIS 

The  natural  frequencies  of  the  flexible  arm  has  been 
determined  experimentally  and  compared  with  analytical 
predictions  obtained  by  using  the  GIFTS  finite  element 
analysis  program.  The  experimental  and  analytical  results 
are  in  good  agreement  except  the  first  mode. 

D.  RECOMMENDATIONS 

1.  The  analytical  model  of  the  experimental  set-up  needs 
to  be  improved  by  full  representations  of  flexible  modes, 
including  control  laws,  and  considering  nonlinear  effects. 
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2.  Future  modal  analysis  should  incorporate  a  filter  to 
reduce  the  effects  of  noise.  Force  windows  and  exponential 
windows  are  available  in  the  data  acquisition  system  for 
such  a  purpose. 

3.  Modal  analysis  using  the  VISTA  data  acquisition 
package  and  IDEAS  modal  analysis  software  should  be 
completed  on  the  structure  once  it  is  connected  together  and 
mounted  on  the  granite  table  on  airpads.  Natural 
frequencies  and  mode  shapes  should  be  then  determined. 
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